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that  processing  characteristics  such  as  fiber  breakage,  overlap  and 
orientation  strongly  Influence  the  fatigue  endurance  of  the  products. 

Both  fatigue  life  and  tensile  strength  can  be  characterized  by  extreme 
value  statistics  and  for  a specific  population  there  exists  a unique  rela- 
tion between  the  cumulative  distribution  of  tensile  strengths  and  the  cumu- 
lative distribution  of  fatigue  lifetimes.  It  was  found  that  for  sifcpopula- 
tlons  derived  from  the  same  main  population,  an  applied  stress -breaking 
time  superposition  procedure  (l.e.  a_av  - log  tp)  can  be  applied,  with  the 

shift  factor  being  a /<o> dependent.  Using  the  master  curve,  one  can  gen- 
erate fatigue  life  distributions  for  different  values  of  0max/<0>  within  a 

reasonable  degree  of  accuracy.  The  environmental  shift  seems  to  be  totally 
expressed  by  the  effect  of  the  environment  on  the  breaking  strength  distri- 
bution. 

Damage  propagation  rate  measurements  Indicate  that  under  a constant 
set  of  upper  and  lower  stress  limits  the  rate  of  damage  Is  constant.  The 
rate  of  propagation  was  found  to  correlate  best  with  a contol nation  of  the 
Initial  K|nax  and  aK  through  the  form 
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Fatigue  life  prediction  through  the  use  of  breaking  strength  distribu- 
tion and  damage  propagation  laws  was  tried.  The  model  proposed  was  able  to 
generate  a reasonable  distribution  curve  at  a single  set  of  experimental 
conditions.  The  damage  propagation  rate  law  had  to  be  empirically  modified 
however.  In  order  to  predict  the  distribution  under  different  sets  of 
experimental  conditions. 

Addition  of  polyaramld  fibers  to  the  graphite  composites  reduced  the 
sensitivity  of  crack  propagation  to  both  the  stress  amplitude  and  the 
maximum  stress.  Hybrid  composites  were  stronger  and  tougher  than  a 
graphite  reinforced  composite  with  an  equivalent  amount  of  fiber  reinforce- 
ment. 
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ABSTRACT 

The  fatigue  failure  of  graphite  fiber  reinforced  nylon  is  Influenced  by  both 

the  technique  of  fabrication  of  the  composite  and  environmental  conditions. 

Compression  molded  materials  exhibited  an  isothermal  brittle  fatigue  failure 

through  initiation  and  propagation  of  cracks.  Injection  molded  materials  were 

non-linear,  exhibiting  ductile  flow  prior  to  ultimate  failure.  It  is  postulated 

that  localized  flow  at  a crack  tip  caused  thermal  melting,  leading  to  a 

thermally  Induced  fatigue  failure.  This  study  Indicated  that  processing 

characteristics  such  as  fiber  breakage,  overlap  and  orientation  strongly  Influence 

the  fatigue  endurance  of  the  products. 

Both  fatigue  life  and  tensile  strength  can  be  characterized  by  extreme 

value  statistics  and  for  a specific  population  there  exists  a unique  relation 

between  the  cumulative  distribution  of  tensile  strengths  and  the  cumulative 

distribution  of  fatigue  lifetimes.  It  was  found  that  for  subpopulations  derived 

from  the  same  main  population,  an  applied  stress-breaking  time  superposition 

procedure  (l.e.  amax  - log  tg)  can  be  applied,  with  the  shift  factor  being 

a_  /^dependent.  Using  the  master  curve,  one  can  generate  fatigue  life 

distributions  for  different  values  of  om.v/<o>  within  a reasonable  degree  of 

max 

accuracy.  The  environmental  shift  seems  to  be  totally  expressed  by  the  effect 
of  the  environment  on  the  breaking  strength  distribution. 

Damage  propagation  rate  measurements  indicate  that  under  a constant  set 
of  upper  and  lower  stress  limits  the  rate  of  damage  Is  constant.  The  rate  of 
propagation  was  found  to  correlate  best  with  a combination  of  the  Initial 
K|nax  and  aK  through  the  form 
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Fatigue  life  prediction  through  the  use  of  breaking  strength  distribution  and 
damage  propagation  laws  was  tried.  The  model  proposed  was  able  to  generate 
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a reasonable  distribution  curve  at  a single  set  of  experimental  conditions. 

The  damage  propagation  rate  law  had  to  be  empirically  modified,  however.  In 
order  to  predict  the  distribution  under  different  sets  of  experimental 
conditions. 

Addition  of  polyaramld  fibers  to  the  graphite  composites  reduced  the 
sensitivity  of  crack  propagation  to  both  the  stress  amplitude  and  the  maximum 
stress.  Hybrid  composites  were  stronger  and  tougher  than  a graphite  reinforced 
composite  with  an  equivalent  amount  of  fiber  reinforcement. 
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FOREWORD 


There  has  been  a growing  use  of  short  fiber  reinforced  polymeric  solids 
In  engineering  components  subject  to  fatigue  failure,  necessitating  a more 
quantitative  representation  of  the  reliability  of  such  structures.  The 
objective  of  this  research  Is  to  develop  for  fiber  reinforced  composites  a 
framework  for  Interrelating  the  statistical  variability  of  the  fracture 
process  with  the  mechanical  and  environmental  history  of  the  structure. 

The  work  described  In  this  report  was  performed  at  the  Institute  of 
Materials  Science,  University  of  Connecticut,  under  Contract  No.  DAAG-46-75-C-0064. 

The  research  was  sponsored  by  the  Army  Materials  and  Mechanics  Research  Center, 

Watertown,  Massachusetts,  with  Dr.  Margaret  Roylance  as  Technical  Supervisor. 
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The  reliability  of  fiber  reinforced  materials  in  engineering  com- 
ponents subject  to  fatigue  failure  Is  of  growing  interest  as  these 
materials  become  more  widely  used.  It  Is  established  that  progressive 
fatigue  damage,  slowly  building  up  to  some  critical  value,  is  the 
principal  mechanism  of  failure  In  fiber  reinforced  plastics  (1-3).  The 
fatigue  process  develops  in  three  stages:  nucleation  of  submicroscopic 
flaws,  slow  growth  of  flaws  under  repeated  load  and  finally  rapid  pro- 
pagation of  cracks  co  c-ltastrophic  failure.  Failure  can  occur  by  either 
thermal  melting  (through  viscoelastic  energy  dissipation)  or  brittle 
failure  (4).  In  this  study  we  will  examine  brittle  fatigue  failure 
only. 

Fiber  reinforced  plastics  possess  a distribution  of  submicroscopic 
and  microscopic  flaws  whose  character  is  determined  primarily  by  the 
fabrication  and  handling  history  of  the  material.  Upon  application  of  a 
loading  cycle,  fatigue  cracks  nucleate  and,  from  this  instant,  the 
fatigue  life  is  controlled  by  their  rate  of  growth.  Although  it  is 
likely  that  initiation  processes  are  important  in. unnotched  polymer 
materials,  there  is  little  doubt  that  the  multitude  of  fiber  ends  and 
debonded  interfaces  in  fiber  reinforced  composites  supply  many  sites  for 
crack  initiation  during  the  very  early  application  of  stress  (5,6). 

In  homogeneous  solids  the  critical  variable  for  characterizing  the 
rate  of  crack  growth  is  the  stress  intensity  factor,  K (7).  For  a 
unidlrectlonally  stressed  elastic  body  containing  a crack  of  length 
(2a),  the  stress  intensity  factor  may  be  written  as: 

K » Y o(ir(a  + ry))1/2  (1) 

where  Y is  a geometrical  parameter  which  depends  on  the  specimen  con- 
figuration (8),  o is  the  applied  stress,  (a)  is  the  crack  half-length 
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and  ry  Is  the  "plastic  ilone"  size.  The  plastic  zone  around  the  tip  of  a 
crack  In  a polymeric  material  Is  the  yielded  or  crazed  material  around 
the  Imbedded  defect  whose  dimensions  depend  on  the  crack  size  and  material 

properties.  Dugdale  (9)  has  shown  that  at  stresses  well  below  the  yield 

JZ  2 

point,  r^  * /Oy ),  where  oy  is  the  yield  stress. 

In  viscoelastic  polymers,  a crazed  zone  develops  around  a crack  tip 
with  time.  The  crazed  zone  can  grow  at  constant  crack  size  until,  at  a 
critical  time,  the  crack  grows  through  the  crazes.  A number  of  theories 
for  delayed  crack  growth  In  viscoelastic  solids  have  been  presented  (10- 
14)  which  express  both  the  rate  of  crazing  and  delayed  crack  growth  as  a 
function  of  stress  intensity  factor  and  time.  A recent  study  of  craze 
and  crack  growth  In  polycarbonate  (14)  Indicates  that  both  processes  can 
be  characterized  by  the  same  time  dependent  function. 

Whether  a stress  Intensity  factor  (developed  for  homogeneous  solids) 
can  be  used  to  correlate  crack  growth  In  a heterogeneous,  viscoelastic, 
fiber  reinforced  composite  Is  certainly  open  to  question.  The  stress 
field  Is  complicated  by  the  presence  of  multiple  phases  and  the  specifi- 
cation of  a fracture  criterion  Is  complicated  by  factors  such  as  fiber 
bridging  at  the  crack  tip.  The  use  of  a single  value  of  critical  stress 
intensity  factor  to  mark  the  onset  of  catastrophic  failure  requires  one 
to  look  upon  the  heterogeneous  body  as  quasi  homogeneous  with  respect  to 


the  fracture  process.  Available  literature  suggests  that  the  breaking 
stress  of  Isotropic  short  fiber  reinforced  materials  Is  Inversely  pro- 
portional to  the  square  root  of  crack  size  as  expected  by  equation  (1) 
(15-17)  and  that  an  "effective"  stress  Intensity  factor  might  be  a 
useful  material  parameter. 

In  the  first  part  of  this  study,  a stochastic  model  relating  the 
cumulative  distribution  of  fatigue  lifetimes  directly  to  the  cumulative 
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distribution  of  breaking  stress  Is  developed  which  by-passes  the  need  of 
defining  a stress  Intensity  factor  or  measuring  critical  flaw  size.  In 
the  second  part  of  the  study  the  rate  of  crack  growth  from  a manufactured 
notch  Is  measured  as  a function  of  "effective"  stress  Intensity  factor 
and  the  resulting  equations  are  used  to  calculate  directly  the  fatigue 
lifetime. 

II.  STATISTICS  OF  BRITTLE  FRACTURE 


Brittle  fracture  generally  originates  In  highly  localized  regions 
where  stress  Is  concentrated  at  submlcroscoplc  flaws.  Fracture  of  a 
specimen  at  a given  level  of  stress  depends  on  the  statistical  expecta- 
tions of  encountering  a critical  (or  weakest)  flaw  having  the  specified 
strength.  If  flaws  are  distributed  at  random  throughout  a population, 
the  distribution  of  strengths  Is  statistically  related  to  the  distri- 
bution of  largest  flaws  appearing  throughout  the  population.  This  model 
Is  Identical  to  the  generalized  asymptotic  theory  of  extreme  values 
(18).  For  samples  consisting  of  n elemental  volumes  having  a probability 
density  function  of  elemental  strength  f(o)  (I.e.  f(o)do  = fraction  of 
elements  having  a strength  between  a and  (o  + do))  and  cumulative  dis- 
tribution of  elemental  strengths  F(o)  = /°f(o)do,  the  distribution  of 
smallest  values  of  o (I.e.  the  fraction  of  specimens  with  a strength 
between  o and  (o  + do))  Is  given  by: 

9„(<0  • "f(o)0  - Ffo))"-1  (2) 

and  the  cumulative  distribution  of  sample  strengths  Is  given  by: 

6n(o)  - /J  g(o)do  - 1 - (1  - F(o))n  OJ 

This  Is  the  probability  of  encountering  a specimen  with  a strength  equal 
to  or  less  than  o.  For  the  case  of  fracture,  the  elemental  probability 
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function  f(o)  Is  unknown  since  the  random  variable  o cannot  be  measured 
directly  except  for  the  extreme  value  Itself  (l.e.  the  fracture  stress 
of  the  material). 

A model  for  the  function  f(a)  has  been  developed  by  Welbull  based 

i 

on  experimental  data  on  the  cumulative  distribution  of  sample  strengths 

G (a)  (19).  A two  parameter  version  of  the  Weibull  function  Is: 

f(o)  - £ (£  exp  - <£  )«  (4) 

0 0 0 

which  results  In  a cumulative  distribution  of  sample  strengths  given  by: 

1 - P(o)  « G(o)  = 1 - exp  - n(f  )“  (5) 

o 

where  aQ  and  a are  constants  and  n Is  the  number  of  elemental  volumes  In 
the  material.  (For  a population  with  specimens  of  fixed  volume,  n Is 
constant  and  may  be  combined  with  the  constant  o0.).  This  type  of 
equation  Is  widely  used  to  characterize  the  variability  of  strength  In 
fiber  reinforced  composites  (20-25). 

The  results  of  fatigue  tests  can  also  be  characterized  by  the 
statistics  of  extreme  values.  The  specimens  that  fail  under  a given 
stress  amplitude  at  various  numbers  N of  load  cycles  may  be  considered 
as  forming  a group  of  the  weakest  specimens  out  of  a large  population  of 
samples.  The  distribution  of  N among  this  subpopulation  may  therefore 
be  described  by  a Welbull  function  of  the  form  (26): 

1 - P(N)S  * G(N)S  - 1 - exp  - (Jj-)°f  (6) 

where  G(N)  Is  the  cumulative  probability  of  encountering  a specimen  with 
a life  of  N cycles  or  less  under  a stress  of  amplitude  S.  The  parameters 
Nq  and  <>f  are  constants. 

Equations  5 and  6 are  only  two  of  a variety  of  cumulative  distribution 
functions  that  might  be  used  to  describe  the  strength  and  fatigue  life 
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distributions  of  fiber  reinforced  materials.  More  complex  forms  must  be 
used  when  factors  such  as  a lower  limiting  strength,  endurance  limits, 
blmodal  elemental  probability  distributions  or  Induction  times  for  crack 
nucleatlon  must  be  considered.  For  the  specific  materials  studied 
herein,  however,  the  above  mentioned  equations  appear  to  be  adequate. 

Regardless  of  the  exact  functional  forms  for  the  cumulative  dis- 
tribution equations.  It  Is  likely  that  if  both  breaking  strengths  and 
cycles  to  failure  can  be  described  as  Poisson  processes,  one  should  be 
able  to  find  a functional  relationship  between  the  two  quantities  at 
equi probability  points  in  the  respective  distributions,  l.e.: 

(N)$  - f(ob)  » P(H)S  - P(ob)  (7) 

This  idea  will  be  used  to  develop  a single  "master-curve"  of  fatigue- 
life  versus  breaking  strength  for  a given  population  of  material.  The 
master  curve  will  then  be  used  to  estimate  changes  In  fatigue  lifetime 
caused  by  environmental  degradation. 

III.  A MODEL  FOR  FATIGUE  CRACK  PROPAGATION 

It  Is  generally  agreed  that  failure  Initiates  from  critical  flaws 
In  the  material  structure.  Thus,  If  one  knows  the  Initial  flaw  size  and 
the  rate  of  flaw  growth  under  a given  loading  history,  one  should  be 
able  to  estimate  the  tlme-to-break.  There  have  been  a variety  of  crack 
propagation  laws  suggested  In  the  literature  (27,28).  An  appropriate 
function  for  crack  propagation  in  fiber  reinforced  plastics  Incorporates 
the  effects  of  a maximum  stress  Intensity  factor  1^  , the  amplitude  of 
stress  (In  the  case  of  fatigue),  a lower  threshold  below  which  no  crack 
propagation  occurs  (Ktb)  and  an  upper  limiting  value  above  which  crack 
propagation  Is  Instantaneous  (K.)  (29).  A model  was  chosen  that  ms  a 

v« 
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combination  of  one  used  by  Erdogan  (30)  for  crack  propagation  In  fatigue 
and  by  Kitagawa  and  Motomura  for  crack  growth  In  viscoelastic  solids 
04): 


where  a Is  the  flaw  grow  rate,  6 Is  a rate  constant  dependent  on  tempera- 
ture, frequency,  environment  etc.,  a Is  a constant  and  ERD  is  the  variable 
defined  by  the  following  equation: 

ERD  ■ ^ 0 - R)t  (9) 

The  variable  is  the  stress  Intensity  factor  at  the  maximum  stress 
In  the  test,  R Is  the  ratio  of  maximum  to  minimum  stress  (or  stress 
Intensity)  during  fatigue  and  y Is  a constant.  The  lower  threshold  is 
characterized  by  ERDT  * IC^  (l-R)^and  the  upper  threshold  Is  given  by 
ERDC  = KgO-R)*. 

Although  It  Is  felt  that  this  Is  a realistic  model  for  crack  growth 
In  the  systems  that  have  been  studied.  It  Is  recognized  that  a five- 
parameter  model  Is  not  required  for  adequate  fitting  of  most  experimental 
data  that  Is  presently  available.  An  equation  of  the  form  (30): 

* ' o-MTir  - *<&  ^ <10> 

although  not  containing  an  upper  and  a lower  bound,  fits  most  data  quite 
satisfactorily  (29). 

Very  little  data  on  composites  are  available,  but  models  similar  to 
the  above  mentioned  are  used  for  correlation  (31-35).  In  orthotropic 
laminates,  resistance  to  crack  propagation  Is  a function  of  orientation 
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and  there  are  a variety  of  energy  dissipating  mechanisms  Including 
fiber/matrix  debonding  to  blunt  crack  propagation  (36-38).  Studies  on 
randomly  dlspersed-short  fiber  reinforced  materials  (15,17)  Indicate 
that  an  equation  of  the  form: 

a = B(aK)°,  (11) 

where  AK  Is  the  amplitude  of  the  stress  and  B and  a are  constants.  Is 
satisfactory  for  fitting  the  experimental  data. 

IV.  FATIGUE  LIFE  PREDICTION 

Several  approaches  to  fatigue  life  predictions  in  composites  have 
been  tried  (39-41).  Broutman  and  Sahu  (39)  used  a modified  version  of 
Miner's  law  to  describe  lifetime  distributions  for  a population  of  E- 
glass  reinforced  epoxy.  For  a similar  system,  Agarwal  and  Dally  (40) 
found  that  fatigue  life  could  be  correlated  to  a normalized  stress 
(Aa)/ou: 

— * a - B log  N (12) 

°u 

where  Ao  is  stress  amplitude  Is  ultimate  strength  and  N Is  cycles  to 
break.  Salee  (29)  and  Hahn  (24)  have  found  unique  relationships  between 
breaking  strength  distributions  and  fatigue  life  distributions.  A more 
general  approach  to  fatigue  life  prediction  has  been  Introduced  by 
Halpln  et.  al.  (41).  They  use  a Helbull  distribution  function  and  a 
power  law  growth  model  to  Obtain  a cumulative  distribution  equation  of 
the  form: 

ob<t)  z(r-,)  . f 

P(oh(t)>o)  ■ exp  - [(-I — ) + (£-)] 

D at 


(13) 
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where  r,  o,  t and  are  constants  and  ab(t)  Is  residual  strength  at 
time  t.  Equation  (13)  can  be  used  to  calculate  the  cumulative  strength 
distribution  at  a specified  time  t **  tQ  or  to  calculate  a cumulative 
distribution  of  breaking  times,  P(tg)  at  a specific  residual  strength 

°b(t)  ■ “max' 

V.  FATIGUE  BEHAVIOR  OF  SHORT  GRAPHITE  FIBER  REINFORCED  NYLON  66  COMPOSITES 
Cumulative  distribution  functions  for  both  tensile  strength  and 
tlme-to-break  In  fatigue  were  generated  for  several  populations  of 
graphite  fiber  reinforced  nylon  66.  In  addition,  crack  propagation  rate 
data  were  obtained.  It  will  be  shown  that  for  a population  of  composite 
materials  there  Is  a unique  relationship  between  tensile  strength  and 
fatigue  life  and  that  a-prlorl  prediction  of  fatigue  lifetime  from  crack 
propagation  rate  data  Is  extremely  difficult  because  of  the  high  sen- 
sitivity of  propagation  rate  to  stress. 

1 . Experimental  Procedures 
Materials 

The  experimental  work  was  performed  on  short  graphite  fiber- 
reinforced  nylon  66  materials.  The  nylon  66  was  a 60-80  mesh  powder  of 
DuPont's  2yte1  101.  The  average  modulll  of  compression  molded  specimens 
of  the  nylon  were  1.8  x 10^  psl  and  7 x 10*  psl  at  55  and  100  percent 
Relative  humidity  respectively.  The  graphite  fibers  were  Hercules  type 
AS,  with  6 nm.  average  length,  8.4  x 10  mm.  average  diameter,  420,000 
psl  average  strength,  32  x 10®  psl  average  modulus  of  elasticity  and  1.8 
grams/cc  (0.065  lb/1n^)  density.  Two  batches  of  fibers,  obtained  one 
year  apart,  were  used.  The  two  batches  were  somewhat  different  In  pro- 
perties and  will  heretofore  be  referred  to  as  fiber  batch  1 and  fiber 

batch  2.  The  compression  molding  fabrication  technique  has  been  reported 

* 

previously  (23). 


1 


' 1}  The  composite* samples  contained  10%  by  volume  of  fiber  in  a random-in-a- 

I t plane  orientation.  Standard  6.5  Inch  tensile  bars  were  machined  from  1/8  Inch 

t thick  composite  plates.  Cantilever  beams  for  use  In  the  crack  propagation 

experiments  were  6 In.  by  1/2  In.  by  1/8  in.  Precisely  placed  channels  were  cut 
in  the  center  of  each  half-inch  face  parallel  to  the  length  forming  a path  for 
crack  propagation.  A sharp  crack  1.0  In.  in  length  was  cut  into  one  side  of  the 
beam  between  the  loading  points.  In  addition  to  the  compression  molded  specimens 
prepared  In  the  laboratory,  a series  of  injection  molded  tensile  bars  containing 
13.7  volume  per  cent  fibers  were  obtained  from  LNP  Corporation. 

The  test  specimens  were  annealed  for  24  hours  In  a silicone  oil  bath  at 
185°C.  They  were  then  conditioned  for  24  hours  In  a boiling  48  per  cent  aqueous 
solution  of  potassium  acetate.  During  this  process  the  composite  absorbs  2.8 
per  cent  water  which  is  equal  to  the  equilibrium  adsorption  in  air  at  50  per  cent 
relative  humidity. 

Mechanical  Testing  Procedure 

Tensile  stress-strain  measurements  were  made  on  an  Instron  Model  1230  ■ 

* 

hydraulically  actuated  testing  machine  at  a crosshead  speed  of  0.2  In.  per  minute 

: 

and  room  temperature.  The  average  cross-sectional  area  In  the  gauge  length  and  \ 

\ 

the  actual  cross-section  at  the  break  were  measured  and  the  elastic  modulus,  1 

| 

breaking  strength,  proportional  limit  and  elongation  at  break  were  determined. 

Cumulative  distribution  curves  for  tensile  fatigue  at  10  cycles  per  second 
l were  also  obtained.  The  maximum  and  minimum  loads  for  the  fatigue  test  were 

r 

i calculated  and  the  machine  set  In  the  load  control  mode.  In  this  mode  the  load 

limits  are  maintained  regardless  of  the  creep  In  the  sample.  In  load  control, 
the  mean  load  and  command  are  first  set  to  zero.  The  function  generator  Is  set 
to  give  a sine  wave  at  10  Hertz  and  then  the  mean  load  Is  Increased  to  the  required 
amount.  The  command  dials  are  Increased  until  the  maximum  and  minimum  load  limits 
* are  obtained.  This  process  usually  took  of  the  order  of  10  seconds  or  roughly  100 

A'  cycles.  At  this  point  the  cycle  counter  Is  started  and  the  test  begins. 
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Flaw  Propagation  Rate  Measurements 

In  the  graphite  fiber  reinforced  materials,  opacity,  crack  bridging 
by  broken  fiber  ends  and  Irregularly  shaped  fracture  surfaces  made  it 
Impossible  to  follow  the  crack  tip  visually.  The  propagation  rate  was 
measured  Indirectly  by  measuring  crack  opening  displacement.  The  pre- 
notched cantilever  beams  were  precall  bra ted  by  machining  cracks  of  a 
given  size  and  measuring  beam  compliance  and  crack  opening  displacement 
as  a function  of  load.  An  effective  stress  Intensity  factor  was  defined 
by  the  usual  definition  In  terms  of  compliance  of  a double  cantilever: 

„2  _ 3P2  [4(a  + 0.6h)2  . 1-v  1 

W"  L i,r  + TTJ  (1/ 

where  P Is  the  applied  force  at  the  loading  pins,  W Is  the  width  of  the 
notched  section,  b Is  beam  width,  h Is  the  beam  height,  v Is  Poisson's 
Ratio  and  a Is  the  flaw  length. 

2.  Experimental  Results  and  Analysis  of  Data 
Data  are  reported  on  four  populations  of  graphite  fiber  reinforced 
nylons  which  will  be  given  the  following  designations 

TABLE  1.  Experimental  Populations 


Designation Treatment 


1 A \ RH  55 

Nylon  66,  Fiber  Batch  1,  compression  molded, 
annealed  1/2  hour,  conditioned  at  55  percent 
relative  humidity 

1 A 1 RH  100 

Nylon  66,  Fiber  Batch  1,  compression  molded, 
annealed  1 hour,  conditioned  at  100  per  cent 
relative  humidity 

2 A 24  R.H.  55 

Nylon  66,  Fiber  Batch  2,  compression  molded, 
annealed  24  hours,  conditioned  at  55  per  cent 
relative  humidity 

LNP  Commercially  available  graphite  reinforced  nylon  66, 

annealed  24  hours  and  conditioned  at  55  per  cent 
relative  humidity 


yasrvf: 
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The  tensile  stress-strain  curves  were  linear  up  to  the  breaking  point 
except  for  37  percent  of  the  specimens  In  the  1 A 1 RH  100  population 
which  showed  some  yielding  at  9500  + 300  psl.  , The  Initial  moduli  of 
elasticity  are  presented  In  Table  2: 


TABLE  2. 


Moduli  of  Elasticity 


Population 


1 A j RH  55 
1 A 1 RH  100 


Modulus  of  Elasticity 




950.000  + 100,000 

795.000  + 145,000 


2 A 24  RH  55 


1,079,000  + 265,000 


The  lower  modulus  In  the  1 A 1 RH  100  specimens  relative  to  that  of  the 
1 A 1/2  RH  55  Is  caused  by  the  lower  modulus  of  water  saturated  nylon  66 
(7  x 104  psl)  compared  to  nylon  66  at  55  percent  humidity  (1.8  x 105  psl). 
The  slightly  higher  modulus  of  the  2 A 24  RH  55  population  Is  caused  by 
a somewhat  higher  degree  of  crystallinity  In  the  nylon  annealed  for  24 
hours.  The  variability  In  all  cases  Is  associated  primarily  with  composi- 
tional variations  from  piece  to  piece.  A one  percent  variation  In  the 
fiber  composition  (reasonably  expected  with  the  layup  procedure  employed) 
can  lead  to  a + 120,000  psl  variation  In  the  modulus.  The  theoretical 

g 

moduli,  calculated  using  the  Pagano-Tsal  equation  (42),  are  1.29  x 10 
psl  and  1.01  x 106  psl  for  the  populations  conditioned  at  55  percent  RH 
and  100  percent  R.H.  respectively.  Fiber  efficiencies  are  thus  of  the 
order  of  75  to  85  percent.  Indicating  good  stress  transfer  from  matrix 
to  fiber.  The  higher  modulus  of  the  LNP  fibers  Is  due  to  the  higher  fiber 
loading  (13.7  per  cent  rather  than  10  per  cent). 
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The  cumulative  distributions  of  tensile  strengths  are  shown  in  Figures  1 
and  2.  All  populations  conform  to  a simple  Weibull  function  with  coefficients 
of  determination  greater  than  98  per  cent.  The  Weibull  parameters,  the  average 
breaking  strengths  and  the  coefficients  of  variation  are  summarized  in  Table  3. 
TABLE  3 Weibull  Parameters  for  Tensile  Strength  Distribution  Function 


Population 

a 

a (psi) 

<o>  psi 

Coeff.  of  Var. 

1A  1/2  RH  55 

9.9 

13,144 

12,409 

0.112 

1 A 1 RH  100 

11.7 

11,446 

10,985 

0.078 

2 A 24  RH  55 

14.8 

16,768 

16,190 

0.078 

LNP 

23.6 

20,097 

19,688 

0.047 

If  we  estimate 

that  the 

theoretical 

strength  of  a 

random-in-a-plane  reinforced 

composite  Is  roughly  equal  to  one-third  of  the  value  predicted  by  the  rule-of- 
mlxtures,  the  predicted  average  strength  should  be  of  the  order  of  16,000 
to  17,000  psl  (based  on  a fiber  strength  of  420,000  psi  and  a nylon  strength 
of  8,300  psl  at  55%  RH  and  6,300  psl  at  100%  RH).  This  leads  to  fiber  utilization 
efficiencies  of  70  to  90  per  cent,  which  are  In  agreement  with  the  values 
calculated  from  the  modulus.  The  Weibull  parameters  for  the  2 A 24  RH  55 
population  were  higher  than  for  the  other  compression  molded  populations.  This 
group  of  specimens  was  prepared  using  a different  batch  of  fibers.  Although 
the  two  were  supposedly  Identical,  It  appears  that  the  fibers  In  the  second  batch 
are  somewhat  stronger  and  have  less  variation  In  their  strength.  (In  retrospect 
It  was  recalled  that  the  handleablllty  of  the  two  fiber  batches  was  also  somewhat 
different).  The  coefficient  of  variation  for  the  LNP  injection  molded  materials 
Is  much  smaller.  Indicating  a more  reproducible  process  for  preparing  the  specimens. 
The  conditions  for  the  cyclic  fatigue  tests  are  summarized  In  Table  4: 


T 
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TABLE  4 Stress  Amplitudes  in  Fatigue  Experiments 


Population Conditions 


wsi 

°mean  psi 

R =(0min/°max) 

amax/<a> 

1 A 1 RH  100-1 

9061 

5458 

0.20 

0.83 

1 A 1 RH  100-2 

9279 

5458 

0.18 

0.85 

1 A 1 RH  100-3 

9824 

5458 

0.11 

0.90 

1 A 1/2  RH  55 

11250 

6250 

0.11 

0.90 

2 A 24  RH  55 

13566 

8172 

0.20 

0.83 

LNP 

13788 

7873 

0.14 

0.70 

Because  the  maximum  stress  in  the  tests  were  close  to  the  average 
strengths  of  the  populations,  a fraction  of  the  weaker  specimens  broke  on  the 
first  half-cycle  (represented  by  G(o)  below).  The  1 A 1 RH  100-1  population 
also  possessed  a considerable  number  of  specimens  surviving  over  10^  cycles 
( 1 05  seconds),  indicating  that  these  composites  have  an  endurance  limit  that  is 
at  a relatively  high  stress  level.  The  Wei  bull  function  used  to  describe  the  data 
was  modified  to  account  for  the  above  mentioned  facts: 

P(tB)  = (l-G(o))  exp  - (tB-t0/t)af  + P(«)  (15) 

The  parameters,  determined  by  a linear  regression  analysis,  are  given  in 
Table  5,  while  plots  of  the  data  are  shown  in  Figures  3 and  4. 

TABLE  5 Wei  bull  Parameters  for  Fatigue  Life  Distribution  Function 


Population 

af 

t (seconds) 

G(o) 

PH 

t„  seconds 
0 

1 A 1/2  RH  55 

0.68 

4661 

0.103 

0.000 

0 

1 A 1 RH  100-1 

0.44 

5852 

0.148 

0.148 

0 

1 A 1 RH  100-2 

0.44 

6576 

0.000 

0.000 

0 

1 A 1 RH  100-3 

0.82 

520 

0.273 

0.000 

0 

2 A 24  RH  55 

0.59 

6034 

0.093 

0.000 

0 

LHP 

0.53 

828 

0.000 

0.000 

298 

i 


4 


i 


I 


- 14  - 


Coefficients  of  determination  are  greater  than  95  per  cent  for  all  cases. 

Crack  propagation  rate  measurements  were  carried  out  on  specimens  from  the 
three  compression  molded  populations.  The  interaction  between  the  stress 
intensity  range,  aK,  and  the  absolute  value  of  the  maximum  stress  intensity,  1^ 
was  determined  by  changing  both  values  during  the  course  of  an  experiment. 

The  most  important  characteristic  of  the  flaw  propagation  is  that 
the  rate  depends  only  on  the  initial  value  of  the  stress  intensity 
parameter.  At  a fixed  value  of  a and  R,  the  flaw  propagation  rate  is 
constant  for  almost  the  whole  life  of  the  specimen.  Only  in  the  last 
few  cycles  of  life  does  the  flaw  growth  rate  accelerate.  Thus,  the 
mechanism  of  failure  is  as  follows:  Stress  concentration  around  the 
flaw  tip  causes  the  region  around  the  tip  to  progressively  decrease  in 
strength.  This  zone  grows  linearly  in  time,  but  the  actual  crack  tip 
does  not  move.  When  the  material  has  been  weakened  enough,  crack  pro- 
pagation occurs  rapidly  and  failure  ensues.  If  the  load  is  increased 
prior  to  catastrophic  failure,  the  crack  opens  and  grows  to  a new 
"equilibrium"  size  and  the  damage  then  occurs  at  a faster  constant  rate. 

This  Is  similar  to  the  mechanism  of  craze  growth  In  amorphous  plastics 
under  stress  corrosion  environments. 

The  Erdogan  equation  (Eqn.  10)  provided  a suitable  function  for 
correlation.  By  letting  the  exponent  y vary  from  0.0  to  1.0,  the  relative 
Influence  of  K|nax  and  aK  could  be  explored.  By  running  a regression 
analysis,  the  values  of  y and  r giving  the  best  fit  of  the  data  can  be 
found.  The  results  for  the  three  populations  are  presented  In  Table  6. 

A typical  plot  of  the  experimental  data  Is  shown  as  Figure  5. 
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TABLE  6. 

Parameters  for  the  Erdogan  Eqn.  (101. 

Population 

B 

V 

r 

Coeff.  of 

(a  In  cm/ml n) 

Det. 

* K„,M00  kg-a»-3/2 

1 A 1 RH  100 

7.68  x 10-42 

0.6 

15.0 

0.69 

1.24  x 10‘3 

1 A 1 RH  55 

4.36  x 10“47 

0.1 

16.4 

0.72 

1.43  x 10‘4 

2 A 24  RH  55 

1.68  x 10’57 

0.4 

19.4 

0.83 

8.98  x 10-8 

The  last  column  of  Table  6 indicates  the  crack  propagation  rates  at 
-3/2 

Snax  = 48®  kg-cm  ' and  R=0.2,  which  are  representative  of  the  values 
used  in  the  fatigue  experiments.  The  1 A 1 RH  1 00  population  is  least 
resistant  to  crack  propagation,  while  the  2 A 24  RH  55  population  is 
most  resistant.  Although  a linear  plot  on  log-log  coordinates  appears 
reasonable  (Figure  5),  the  coefficients  of  determination  are  rather  low. 
A plot  of  the  data  on  an  expanded  scale  for  the  1 A 1 RH  100  population 
is  shown  in  Figure  6 and  suggests  the  source  of  the  low  correlation. 

The  data  appear  to  indicate  both  a lower  and  upper  threshold  for  crack 
propagation  as  suggested  by  equation  (8),  but  the  scatter  of  experi- 
mental results,  associated  with  uncertainties  In  the  measured  values  of 
both  velocity  and  maximum  stress  Intensity  make  it  difficult  to  find  a 
"better"  model.  A grid  search  technique  was  used  to  evaluate  the 
parameters  In  equation  (8).  Results  for  two  systems  are  presented  in 
Table  7. 


TABLE  7.  Parameters  for  Equation  (8)  . 


6 

ERDC 

ERDT 

a In  cm/ml n 

k 9 

' rm-3 n 

- cm  y 

a 

C5e7T 

of  Det. 

1 A 1/2  RH  55 

0.93  x 10'8 

635 

425  0.1 

1.0 

0.70 

1 A 1 RH  100 

4.95  x 10"12 

560 

315  0.6 

1.7 

0.71 

1 A 1 RH  100 

1.42  x 10’16 

521 

354  0.6 

2.1 

0.68 

(based  on  fatigue  data) 
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Coeffl clents  of  determination  were  of  the  same  order  as  before,  in 
spite  of  two  additional  parameters.  (The  third  entry  In  Table  7 will  be 
discussed  later).  When  crack  propagation  rate  is  so  highly  sensitive 
to  stress  intensity  factor,  it  will  be  very  difficult  to  develop  models 
with  a greater  degree  of  certainty  than  those  presented  above. 

3.  Fatigue  Life  Prediction:  Correlation  Between  Tensile  Strength 
and  Fatigue  Life 

The  cumulative  distributions  of  both  static  tensile  strength 
and  fatigue  lifetimes  have  been  represented  by  Welbull  functions.  As 
suggested  by  equation  7,  one  should  be  able  to  find  a functional  relation 
between  these  two  variables  at  equivalent  points  In  their  respective 
distributions.  Using  the  data  In  Tables  3 to  6,  one  can  generate  the 
following  master  curve  for  all  the  data  obtained  on  composites  contain- 
ing fibers  of  batch  1: 

t o U 

log,0  ob  > 0.112  [-£*  (1-R)0  SJ  log10  t„  + 3.915  (16) 

The  result  Is  plotted  In  Figure  7.  The  coefficient  of  determination 
for  the  straight  line  Is  0.99.  Thus,  all  the  data  for  a given  population 
of  composites,  made  by  the  same  fabrication  process  and  using  the  same 
raw  materials,  can  be  reduced  to  a single  master  curve.  The  effects  of 
stress  and  environment  on  the  fatigue  life  can  therefore  be  determined 
from  the  effects  on  the  distribution  of  tensile  strengths. 

The  relation  for  the  composites  prepared  using  fibers  of  batch  2 


log10  ab  * 0,101 


o-i<>0-6]' 


1o9l0  *B  + 4,114 


The  curve  Is  nearly  parallel  to  the  former,  but  displaced  to  lower 
values  on  the  time  scale  (Figure  7).  The  slopes  of  these  plots  are 
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associated  with  the  dispersion  of  the  cumulative  distributions  and  are 
thus  related  to  the  distribution  of  flaws  In  the  two  populations.  Since 
the  distribution  of  flaws  Is  probably  controlled  mainly  by  the  fabrica- 
tion process,  the  two  populations  should  exhibit  a parallel  response. 

Thp  "Intercepts"  are  controlled  by  both  the  scale  and  dispersion  of  the 
distributions  and  are  clearly  different  for  the  two  fiber  batches.  In 
this  case,  the  added  fiber  strength  does  not  translate  Into  a propor- 

i 

tlonal  Increase  in  fatigue  life. 

4.  Fatigue  Life  Prediction  From  Crack  Propagation  Data 

The  wearout  model  (41).  equation  13,  was  tested  as  a tool  for 

fatigue  life  prediction.  The  experimental  data  were  used  in  a grid 

search  technique  and  a linear  regression  procedure  to  determine  the 

constants  r,  t and  (o__v/o)2^r"1 ^ of  equation  (13).  The  results  are 

max 

summarized  In  Table  8. 

TABLE  8.  Parameters  for  Wearout  Model  (Eqn.  13) 


Population 

r 

a 

t(sec) 

<w;)2(r-1) 

Coeff.  of 
Determination 

1 A 1/2  RH  55 

11.5 

3116 

-4100 

0.96 

1 A 1 RH  100-1 

22.8 

919 

-0.10 

0.98 

1 A 1 RH  100-2 

16.5 

1060 

0.16 

0.96 

1 A 1 RH  100-3 

18.8 

102 

-0.99 

0.98 

2 A 24  RH  55 

19.0 

365 

-0.12 

0.98 

Although  the  model  fits  the  data  very  well,  the  quantity  r Is  not  In- 
dependent of  environmental  conditions,  as  assumed,  and  the  "best-fit" 
value  of  (°m«x/o)2^r'^  Is  negative  In  many  cases,  which  Is  not  physically 
correct.  It  Is  quite  possible  that  a three  parameter  optimization  could 
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lead  to  other  "best-fit”  values,  but  It  Is  felt  that  these  problems 
render  the  model  much  less  useful  as  a tool  for  fatigue  life  prediction 
for  these  materials. 

A prediction  of  the  breaking  time  In  fatigue  directly  from  the 
crack  propagation  rate  equations  was  also  attempted.  The  procedure  is 
described  below  and  illustrated  in  Figures. 

Start  with  the  distribution  curve  for  1 A 1/2  RH  55.  Since  o_  is 

max 

11,250  psi,  specimens  with  breaking  strengths  less  than  or  equal  to  this 
value  must  have  zero  breaking  time.  Thus,  from  Table  7 and  equation  8, 
ERD  -*•  635  (l.e.  a -*•  <•>)  at  the  point  on  the  distribution  corresponding  to 
ob  ■ 11,250  psi  (at  P (ob)  = 0.8).  The  value  of  the  parameter  (ERD),  and 
thus  the  crack  propagation  rate  a,  car.  be  calculated  at  any  other  point 
on  the  distribution  since  (ERD)^/(ERD)2  * ab2^°bl'  ^ p(*g)  * P(°t,)  = 
0.5,  for  example,  ERD  * 563  and  a * 0.012  cm/min.  Using  this  point  as  a 
reference  point,  one  may  then  calculate  the  relative  crack  propagation 
velocity  at  any  other  point  on  the  distribution.  Since  the  growth  rates 
are  constant  for  most  of  the  life  of  the  material,  one  may  assume  the 
relative  breaking  times  are  Inversely  proportional  to  the  relative 
velocities  (l.e.  *l/a2  * *82^*81^*  The  solid  curve  In  Figure  8 has  been 
generated  In  this  fashion  and  follows  the  exprlmental  distribution  very 
closely.  In  a similar  manner,  one  may  use  this  reference  point  and  the 
crack  propagation  data  for  the  other  populations  to  generate  theoretical 
distributions  of  fatigue  life  at  other  conditions.  This  was  attempted 
using  the  "best-fit"  values  for  the  crack  propagation  rate  (Table  7). 

The  predicted  distribution  was  too  sharp  (Eqn.  1 of  Figures).  The 
calculations  were  redone  using  a different  set  of  values  In  the  crack 
propagation  equation  (3rd  line  of  Table  7 and  Figure *6)  with  a high 
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degree  of  success.  The  solid  lines  in  Figures  4 and  8 were  generated  In 
this  manner. 

Unfortunately,  we  cannot  find  any  a-priori  justification  for  using 
the  latter  set  of  parameters  to  describe  the  crack  propagation  data  for 
the  1 A 1 RH  100  population,  other  than  the  fact  that  they  are  consistent 
with  the  breaking  time  results.  Because  of  the  high  sensitivity  of 
crack  propagation  to  stress  level  and  the  Inability  to  make  more  precise 
measurements  of  the  damage  rate,  predictions  of  lifetime  directly  from 
the  propagation  rate  data  must  be  recognized  as  highly  unreliable  for  this  system. 

5.  Mechanism  of  Fatigue  Failure 

It  was  observed  that  damage  propagation  rates  remained  constant  for 
a fixed  level  of  upper  and  lower  cyclic  stresses  and  that  as  the  stress 
levels  were  changed,  the  damage  propagation  rate  also  changed.  This 
suggests  that  the  propagation  rate  depends  on  the  stress  Intensities 
operating  at  the  tip  of  the  Initial  crack  and  that  the  damaged  region 
around  the  crack  tip  grows  (or  decreases  in  strength)  In  such  a way  that 
the  effective  magnitude  of  stress  concentration  at  the  crack  tip  does 
not  change.  A similar  behavior  was  reported  for  environmental  craze 
growth  In  PMMA  (43). 

Owen  and  Dukes  (44)  observed  that  In  random  chopped  strand  E-glass 
reinforced  laminates  the  Individual  strands  lying  perpendicular  to  the 
stress  direction  were  the  first  sites  of  damage.  The  damage  took  place 
In  the  form  of  debonding  (l.e.  separation  between  the  fiber  and  the 
matrix)  that  Intensified  progressively,  affecting  fibers  at  smaller 
angles  to  the  applied  load  until  at  some  hlghdr  load  (less  than  the 
ultimate  load)  resin  cracks  occurred  In  the  resin  rich  zone  of  the 
material. 
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In  the  nylon -graphite  composite,  the  early  damage  also  occurs  near 
fibers  lying  perpendicular  to  the  stress  direction.  However,  as  can  be 
seen  In  Figure  9,  Instead  of  fiber  debonding  (which  may  also  occur  to  a 
small  extent),  good  adhesion  to  the  fibers  causes  the  resin  to  draw  In 
the  direction  of  the  stress.  The  high  degree  of  elongation  causes  the 
fibers  In  this  region  to  break  randomly.  Broken  ends  may  provide  nuclei 
for  the  concentric  strlatlons  pictured  In  Figure  10.  This  process  may  be 
accompanied  by  cavitation  of  the  resin  about  the  fiber,  which  Involves 
drawing  of  the  resin  In  a direction  radial  to  the  fiber.  Fibers  which 
break  prior  to  matrix  failure  have  to  be  pulled  out  during  the  crack 
opening  process.  The  length  of  fiber  that  has  to  be  pulled  out  (l.e. 
half  the  critical  fiber  length  lc/2)  Is  estimated  to  be  0.12  to  0.17  mm. 
The  breaking  of  fibers,  the  pull  out  of  fibers  bridging  the  damaged 
matrix  and  the  drawing  and  cavitation  of  the  resin  about  the  fibers  are 
the  major  energy  dissipating  processes  occurring  In  the  region  around 
the  crack  tip.  The  damage  zone  around  the  crack  tip  grows  linearly  with 
time  while  the  actual  crack  tip  Is  stationary. 

This  Is  phenomenologically  similar  to  the  drawing  and  cavitation  of 
polymer  bridging  the  faces  of  a propagating  craze.  Once  the  damage  zone 
reaches  Its  critical  dimensions,  a fast  crack  propagates  at  a rate  that 
does  not  allow  further  plastic  deformation  to  occur. 


| 
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Populations 

One  can  see  that  the  average  strength  of  the  injection  molded  specimens  is 
higher  than  that  of  the  compression  molded  materials  and  that  the  cumulative 
distribution  for  the  injection  molded  materials  is  narrower  than  that  of  the 
compression  molded  population.  The  lower  coefficient  of  variation  shows  the 
higher  degree  of  reproducibility  of  the  injection  molding  technique . (Figure  2). 

The  cumulative  distribution  of  breaking  times  for  two  populations  tested 
under  approximately  the  same  conditions  of  stress  are  plotted  in  Figure  3. 

Although  the  upper  limits  of  the  stress  amplitude  applied  to  both  populations  are 
almost  the  same,  the  injection  molded  specimens  were  fatigued  at  only  70%  of  their 
average  strength  while  the  compression  molded  samples  were  fatigued  at  83%  of  the 
average  strength.  In  spite  of  this,  the  mean  breaking  time  (measured  at 
P(tg)  = 0.5)  is  680  seconds  for  the  injection  molded  compared  to  2700  seconds  for 

the  compression  molded  materials.  The  compression  molded  specimens  have  an 
average  fatigue  life  clearly  superior  to  the  injection  molded  specimens.  On 
the  other  hand,  the  compression  molded  specimens  show  a much  broader  variation 
in  fatigue  life  and,  in  fact,  the  lower  20%  of  the  distribution  lies  below  that 
of  the  injection  molded  materials.  Thus,  although  the  average  properties  are 
superior,  one  would  have  to  expect  a greater  number  of  early  failures  in  the 
compression  molded  populations.  In  this  case,  "reliability"  is  not  related  to 
"average  lifetime". 

Typical  temperature  measurements  as  a function  of  time  are  shown  in  Flgurell. 

In  all  cases,  the  surface  temperature  of  the  injection  molded  specimens  Increased 
steadily  with  time  until  the  specimen  failed  at  a temperature  in  the  range  of 
70-80°C.  On  the  other  hand,  the  temperature  of  the  compression  molded  specimens 
reached  a steady  state  of  approximately  30-35°C  after  a few  hundred  seconds. 
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This  suggests  that  the  macroscopic  ductility  of  the  Injection  molded 
materials  led  to  a viscous  heating  of  the  material  and  that  the  failure  mechanism 
was  at  least  partially  influenced  by  the  thermal  effect.  The  induction  time 
tQ  = 298  sec.  is  most  likely  the  minimum  time  required  to  heat  the  specimen  to  a 
critical  temperature  under  the  test  conditions.  The  compression  molded  specimens, 
on  the  other  hand,  are  linearly  elastic  to  the  point  of  failure  and  heating  is  due 
only  to  viscoelastic  energy  dissipation.  The  failure,  in  this  case,  is  an 
isothermal,  brittle  fracture. 

Scanning  electron  micrographs  of  fracture  surfaces  from  the  two  populations 
are  shown  in  Figures  12  and  13.  The  surfaces  of  injection  molded  specimens  clearly 
show  thermal  melting  and  a subsequent  dewetting  of  the  fibers  in  the  slow  crack 
region.  Photographs  of  the  fast  crack  region  confirm  that  the  adhesion  was  initially 
present.  The  fracture  surfaces  of  the  compression  molded  materials,  on  the 
other  hand,  show  good  adhesion  at  all  points.  Examination  of  suspensions  of 
recovered  fibers  from  injection  molded  specimens  indicates  considerable  damage  ! 

i 

and  a wide  distribution  of  fiber  lengths  (Figure  14).  Fiber  lengths  vary  between  | 

I 

0.01  and  1.0  mm  with  an  average  length  being  of  the  order  of  0.1  to  0.3  mm.  | 

The  average  aspect  ratio  is  thus  about  18.  The  fibers  also  exhibit  a longitudinal  | 

orientation.  This  was  determined  by  dissolving  away  nylon  at  the  surface  of  the  j 

material  and  examining  the  exposed  fibers  by  dark  field  photography  (Figure  15). 

■ 

This  preferred  orientation  is  undoubtedly  contributing  to  the  apparent  high 

fiber  efficiency.  The  compression  molded  specimens,  on  the  other  hand,  show  random  | 

orientation  and  no  damage  from  processing.  1 
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VI.  FATIGUE  BEHAVIOR  OF  POLYARAMI D/GRAPHITE  FIBER  REINFORCED  NYLON  66  COMPOSITES 

1.  Fabrication  of  Composite  Test  Specimens 

The  experimental  work  was  performed  on  four  types  of  composite. 

1.  A 20  per  cent  by  volume  short  graphite  fiber  reinforced  Nylon  66. 

2.  A 20  per  cent  by  volume  short  polyaramid  fiber  reinforced  Nylon  66. 

3.  A 10  per  cent  short  graphite/10  per  cent  polyaramid  fiber  reinforced 
Nylon  66,  hybrid  (20  per  cent  fiber  by  volume)  with  the  fibers  mixed 
together  to  form  a random  dispersion. 

4.  A 10  per  cent  short  graphite/10  per  cent  polyaramid  fiber  reinforced 
Nylon  66,  hybrid  (20  per  cent  fiber  by  volume)  laminated  with  alternate 
layers  of  the  two  fibers. 

The  nylon  was  a 60-80  mesh  powder  of  DuPont's  Zytel  101,  supplied  by  Liquid 
Nitrogen  Processing  Corporation.  The  graphite  fibers  were  Hercules,  Inc.,  type  AS, 
(see  page  8).  The  polyaramid  fibers  were  DuPont's  Kevlar-49  type  III,  with 
12.7  mm  average  length.  The  average  diameter  was  1.016  x 10  mm,  the  average 
ultimate  tensile  strength  was  400,000  psl,  the  modulus  of  elasticity  was 
19  x 106  psl  and  the  density  was  1.45  g/cc. 

The  composite  samples  contained  20%  by  volume  of  fiber  In  a random  In-a- 
plane  orientation  and  measured  approximately  7.5  x 7.5  x 1/8  In.  Polyaramid 
fibers  were  slurried  In  a tank  of  water  using  a hand  mixer*  The  dilute  mixture 
was  stirred  for  three  minutes  and  then  allowed  to  rest  for  five  minutes  to  permit 
settling  of  the  fibers  In  the  slurry.  A rigid  net  at  the  bottom  of  the  tank  was 
then  pulled  upward  catching  the  floating  fibers  and  building  up  a mat  of  fibers 
with  a random  In-a-plane  dispersion.  A piece  of  filter  paper  was  then  placed 
over  the  fibers  on  the  net  and  the  mat  was  dried  In  place.  The  net  was  then 
separated  from  the  filter  paper,  leaving  a thin  uniform  random-ln-a-plane  veil 
of  fibers.  The  essential  feature  of  this  operation  Is  that  the  original  polyaramid 
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fiber  bundles  of  the  chopped  roving  are  broken  up  so  that  Individual  fibers  are 
liberated  and  are  free  to  act  as  individual  elements  of  reinforcement.  Without 
this  process  of  liberation  the  composite  properties  are  poor  and  the  products 
are  much  less  reproducible.  One  layer  in  the  final  composite  material  was  made 
from  12  such  veils  cut  to  fill  an  area  of  19  x 19  cm.  One  portion  of  nylon  powder 
was  then  sprinkled  evenly  over  the  resulting  layer.  The  procedure  was  repeated, 
placing  one  composite  layer  on  top  of  the  next  until  the  proper  thickness  of 
composite  sheet  was  attained.  In  order  to  make  a hybrid  laminated  sandwich 
structure  the  layers  of  graphite  fiber  and  PRD  fiber  were  alternated.  For  the 
hybrid  mixture,  each  layer  consisted  of  a mixture  of  both  graphite  and  polyaramid. 
The  composite  mats  were  then  compression  molded  using  the  same  technique  reported 
previously  (23). 

The  composites  were  annealed  between  heavy  metal  plates  in  a silicone  oil 
bath  at  185°C  for  24  hours.  In  the  process,  internal  stresses  were  relieved 
and  the  nylon  developed  a fixed  percentage  of  crystallinity. 

Cantilever  beams  for  use  in  the  crack  propagation  experiments  were  3 inches 
by  1 inch  by  1/8  inch.  Precisely  placed  channels  were  cut  in  the  center  of  each 
face  parallel  to  the  length  of  the  specimen,  forming  a path  for  crack  propagation 
with  a width  of  0.01  inches.  A sharp  crack  1.0  inch  in  length  was  cut  into  one 
side  of  the  beam  between  the  loading  points. 

The  standard  tensile  bars  for  measurement  of  the  breaking  strength  and  the 
modulus  of  elasticity  were  used.  They  were  6.5  Inches  long  by  1/2  Inch  wide. 

The  ends  of  the  dog-bone  shape  were  3/4  inch  wide. 

The  test  specimens  were  boiled  for  24  hours  in  a 48%  aqueous  solution  of 
potassium  acetate.  During  this  process  the  composite  absorbed  2.8  per  cent  water 
which  Is  equal  to  the  equilibrium  absorption  In  air  at  55  per  cent  relative 
humidity.  The  equilibrium  moisture  content  can  be  controlled  at  any  other  level 
by  changing  the  salt  concentration.  The  specimens  were  then  stored  at  55  per  cent 
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relative  humidity  and  room  temperature  In  a glove  box  until  they  were  tested. 

Test  specimens  conditioned  at  100  per  cent  relative  humidity  were  boiled  In 
pure  water  for  30  hours  and  then  stored  under  water. 

2.  Crack  Propagation  Rate  Measurements 

The  compliance  method  was  used  to  monitor  crack  length.  The  calibration 
was  carried  out  by  Increasing  the  crack  length  In  small  steps  of  0.02  to  0.035  of 
an  Inch  with  a jeweler's  saw  and  recording  the  crack  opening  versus  the  applied 
load.  The  slope  of  each  curve  Is  the  compliance  for  this  specific  crack  length. 
The  calculated  compliance  was  plotted  versus  the  crack  length  as  In  Figures 
B-l  to  B-4  In  Appendix  B. 

In  order  to  record  the  crack  opening  at  the  points  of  load  application 
during  the  cyclic  fatigue  crack  propagation  experiment,  a sampling  and  hold 
circuit  was  Installed  Into  the  peak  monitoring  section  of  the  Instron  1230 
Stress  Function  Generator.  This  sampling  and  hold  device  enabled  continuous 
monitoring  of  the  upper  or  lower  peaks  of  the  cyclic  strain  on  a Houston 
strip  chart  recorder. 

3.  Data  Analysis:  Calculation  of  Stress  Intensity  Factor  and  Crack 

Propagation  Rate 

Continuous  monitoring  of  the  upper  limit  of  stroke  at  the  point  of  load 
application  was  obtained  by  feeding  the  stroke  output  signal  Into  the  sampling 
and  hold  device.  The  ratio  of  the  upper  limit  of  the  stroke  to  the  peak  of  the 
load  amplitude  could  be  calculated  to  produce  the  Instantaneous  compliance 
C - ^max  (18) 


The  change  of  the  compl lance  with  time  for  a specific  cyclic  loading  condition 
was  calculated  by: 
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(19) 


dC  _ C2  ~ C1 
3t  " tg  - tf 


where  t Is  the  time. 

The  rate  of  crack  propagation  da/dt  can  be  calculated  from  the  compliance 
curve  and  the  experimental  data  from  the  following  equation 
dC  _ dC  v da 

ar  ■ a? x at 

or  (20) 


da  _ dC  . dC 
at  " 3t ' 3a 

where  dC/da  was  calculated  from  the  derivative  of  the  calibration  curve.  The 
derivative  was  obtained  by  fitting  the  data  to  a polynomial  of  the  form: 

C = Aq  + A1  X ai1/2  + A2  x ai1'5  + . . . + A5  x a^*5  (21) 

The  constants  were  obtained  using  a Linear  Regression  Algorithm  and  then 
di fferentiatlng: 

dC/da  = 1/2  x A1  x a^172  + 1.5  x A?  x a^72  + . . . 

(22) 

...  4.5  x A5  x a^*5 

The  strain  energy  release  rate  G can  be  calculated  from: 


G 


1 

W 


,2  dC 
3a 


(23) 


where  W Is  the  width  of  the  crack  (which  Is  narrower  than  the  sample  because 
of  the  machined  grooves)  and  P Is  the  applied  load.  Maximum  and  minimum  values 
of  G,  1e.,  G(nax  and  Gm^n»  can  be  calculated  by  using  the  maximum  and  minimum 
values  of  the  applied  cyclic  load,  P(nax  and  Pm^n. 

The  plane  strain  stress  Intensity  factor  K can  then  be  computed: 

K ■ [G  E / (1  - v2)  ]1/Z  (24) 
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4.  Mechanical  Testing  Results 
(a)  Tensile  Properties 

Tensile  properties  are  reported  for  four  different  composites,  which  will 
be  given  the  following  designations: 

20GR55RH  - Nylon  66  containing  20  per  cent  by  volume  graphite  fibers 
annealed  for  24  hours,  conditioned  at  55  per  cent  relative 
humidity 

PR55RH  - Nylon  66  containing  20  per  cent  by  volume  polyaramid  fibers 
annealed  for  24  hours,  conditioned  at  55  per  cent  relative 
humidity 

HYL55RH  - Nylon  66  containing  20  per  cent  by  volume  graphite  and 

polyaramid  fibers,  (each  fiber  10%),  annealed  for  24  hours, 
conditioned  at  55  per  cent  relative  humidity 
HYL100RH  - Same  composition  as  HYL55RH,  annealed  for  24  hours  and 
conditioned  at  100  per  cent  relative  humidity 
The  theoretical  and  experimental  breaking  strength  and  modulus  of  elasticity 
are  compared. in  Table  (9).  The  Tsal-Halpin  equations  (42)  were  used  to  calculate 
the  theoretical  modulus  from  constituent  properties: 


3/8  E, 


+ 5/8  E, 


where 


1 + 2 (1/D)f  Vf  Ef  + 2(  1/D) f Q-Vf)  En 


>f+  Vf 


Ef(l  + 2Vf)  + 2(1 -Vf ) E„ 


rf)  + 2 + vf 


(1/D)f  is  the  length  to  diameter  ratio  (l.e.,  the  aspect  ratio  of  the 
fibers.)  Em  is  the  modulus  of  the  matrix.  Ef  and  are  the  modulus  and  volume 
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fraction  of  the  fibers,  respectively. 

For  the  hybrid  composites  the  theoretical  modulus  is  assumed  to  be  the 
volume  fraction  average  of  the  two  components.  A schematic  diagram  of  the  hybrid 
composite  is  shown  in  Fig.  16.  In  each  layer  the  fiber  volume  fraction  is  0.2 
while  the  overall  volume  fraction  for  each  type  of  fiber  is  equal  to  0.1. 


E^-j  and  E22  were  calculated  using  the  following  equations: 

1/E22  =TTirGr+  TTT  (28) 

^ com  u22'  br  U22JPRD 

<En>coo, ' <EnVx0-5t  <EnVd x °-5  <29> 


i 


The  properties  of  the  constituents  are: 

Ef  (graphite)  = 32  x 106psi  (1/D)f  (graphite)  = 714 

Ef  (Kevlar)  = 19  x 106  psi  (1/D)f  (Kevlar)  = 1087 

Em  (Nylon)  (55%  RH)  = 1.8  x 105  psi 
Em  (Nylon)  (100%  RH)  = 7.4  x 104  psi 

There  are  no  reliable  theories  for  predicting  the  theoretical  strength  of 
random- in-a-plane  reinforced  composites,  but  if  we  estimate  that  it  should  be 
roughly  equal  to  one-third  of  the  value  predicted  by  the  rule  of  mixture: 

S = 1/3  (S fVf  + Sm  (1-Vf))  (30) 

We  can  get  "theoretical"  results  for  strength. 

The  properties  of  the  constituents  are: 

Sm  (55%  RH)  = 8.3  Kpsl 
Sm  (100%  RH)  * 6.3  Kpsi 
Sf  (graphite)  * 420  Kpsl 
Sf  (PRD)  = 400  Kpsl 

The  measured  composite  efficiency  Is  defined  as  the  ratio  of  the  experimental 
to  theoretical  value  of  either  the  strength  or  modulus. 
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; = efficiency  = — x 100 

Ltheor 

s {31) 

or,  ? = efficiency  = — x 100 

5theor 

These  values  are  tabulated  in  Table  10. 


TABLE  10.  THE  STRENGTH  AND  MODULI  EFFICIENCIES  OF  THE  COMPOSITES 


Material 

Strength  Efficiency 

Modulus  Efficiency 

GR55RH 

60.8% 

50.8% 

PR55RH 

63% 

76% 

HYL55RH 

66.7% 

65.6% 

HYL100RH 

54.5% 

54.8% 

The  composite  efficiency  is  of  the  order  of  50  to  70  per  cent,  which  Is 
typical  for  short  fiber  reinforced  composite  with  fiber  loadings  of  this  order. 

The  composite's  strength  are  ranked  In  the  following  order, 

HYL55RH  > GR55RH  * PR55RH  > HYL100RH 
The  following  order  is  valid  for  the  modulus  as  well. 

(b)  Fatigue  Crack  Propagation 

Cyclic  fatigue  tests  were  run  on  the  same  four  materials  mentioned  before. 

All  samples  were  run  at  6 Hz.  The  damage  propagation  measurements  were  carried 
out  under  humidity  conditions  Identical  to  those  used  for  conditioning  the  samples, 
since  it  was  observed  that  humidity  changes  greatly  Influenced  the  propagation 
rate:;.  Samples  conditioned  at  100  per  cent  relative  humidity  were  tested  while 
over  water.  In  a tank  especially  designed  for  this  purpose. 

Since  the  Interaction  of  both  the  stress  Intensity  range,  aK,  and  the 
absolute  value  of  the  stress  Intensity  factor  as  represented  by  K|nax  was  sought, 
both  values  were  changed  during  the  course  of  the  experiment. 

An  equation  of  the  form  suggested  by  Erdogen  provided  a suitable  functional 
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formula  for  correlation  of  the  data.  This  equation  can  be  written  as: 
a = A Kmaxm  AKn  = A [Kmax  (l-R)n/n+m]mfn  (32) 

Taking  the  logarithms  of  both  sides  we  obtained 

log  a = log  A + m log  + n log  AK  (33) 

A regression  analysis  was  used  to  obtain  the  best  values  of  m,  n and  A.  Results 

for  the  four  materials  tested  are  given  below: 

GR55RH:  a (cm/mln)  = 0.98  x 10"51  [Kmax  (1-R)0,12]  17,87  (34) 

with  a correlation  coefficient  0.924 
PR55RH:  a (cm/min)  = 4.2  x 10-9  [K^  (1-R)0,79]  2,44  (35) 

with  a correlation  coefficient  0.915 
HYL55RH:  a (cm/min)  = 1.31  x 10-20  [K^  (1-R)0,25]  6,5  (36) 

with  a correlation  coefficient  0.804 
HYL100RH:  a (cm/min)  = 8.18  x 10'13  [K^  (1-R)0,47]  2,85  (37) 

with  a correlation  coefficient  0.908 

*Snax  ls  in  un^ts  of  k9  " cm"3^- 

For  a stress*1ntens1ty  factor  value  of  600  kg  - cm’372  and  a value  of  R = 0.5,  j 

which  Is  representative  of  the  values  used  In  the  fatigue  crack  propagation  5 

| 

experiments,  the  crack  propagation  rates  are  calculated:  t 


GR55RH 

9.8  x 10"3  (cm/min) 

PR55RH 

6.6  x 10“3  (cm/min) 

HYL55RH 

4.9  x 10"3  (cm/min) 

HYL100RH 

2.7  x 10"5  (cm/min) 

which  rank  the  resistance  of  the  composites  to  crack  propagation  in  the 
following  order. 

HYL100RH  > HYL55RH  « PR55RH  > GR55RH 
The  propagation  rates  are  also  plotted  as  a function  of  (1-R)  at 
Kni#x  * 645  kg  - cm"372  and  as  a function  of  KJMX  at  R ■ 0.7  in  Figures  17  and 


18  respectively. 

One  may  conclude  the  following: 

1.  The  sensitivity  of  crack  propagation  rate  to  the  amplitude  of  stress 
intensity  (as  represented  by  1-R)  increases  as  the  material  becomes  more  ductile. 

2.  The  resistance  of  PRD  fiber-reinforced  nylon  to  crack  propagation  is  trans- 
lated to  the  hybrid  composite.  The  material  is  stronger  and  tougher  than  the 
graphite  reinforced  composite  with  an  equivalent  amount  of  fiber  and  Is  about 
as  stiff  as  the  graphite  reinforced  material. 

3.  The  sensitivity  of  crack  propagation  to  1^  increases  as  the  composite 
becomes  more  brittle. 

3.  Fractography  Studies 

The  Cambridge  Stereoscan  Scanning  Electron  Microscope  was  used  for  fracture 
studies.  The  fracture  surfaces  produced  in  the  two  different  experiments  (i.e. 
tensile  tests  and  fatigue  crack  propagation  tests)  were  examined.  Samples  were 
first  vacuum  coated  with  a thin  layer  of  carbon  followed  by  a thin  layer  of 
gold,  In  order  to  prevent  static  charge  build-up. 

Specimens  fractured  In  constant  strain-rate  tension  exhibited  a relatively 
brittle  fracture  as  shown  In  Figures  19  to  21.  Fibers  extending  out  of  the 
fracture  surface  are  relatively  short  for  graphite  but  relatively  long  for  PRD. 

The  pulled  out  fibers  are  mainly  those  which  are  not  oriented  parallel  to  the  load 
direction.  The  pullout  length  to  diameter  ratios  (calculated  from  pulled  fibers 
that  are  oriented  parallel  to  the  load  direction)  are  6 and  IS  for  graphite  and 
PRD  respectively. 

If  one  assumes  that  the  critical  aspect  ratio  Is  of  the  order  of  twice  this 
value,  one  can  estimate  the  shear  strength  at  the  Interface  using  the  shear  lag 
analysis. 
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T = 2I 17*777 


For  graphite  fiber 


_ 420,000  _ ,,  r 
t - ' 2n71~"  ” 17.5  Kpsl 


For  PRD  fiber 


’ • Tfjsf  ' 6-67  "p51 


The  shear  strength  between  graphite  and  nylon  is  of  the  order  of  the  shear 
strength  of  the  nylon  matrix.  This  indicates  that  the  adhesion  of  the  nylon  to 
the  graphite  fibers  is  good.  On  the  other  hand,  the  adhesion  of  the  nylon  to  the 
PRD  fibers  is  poor.  Eagles  (45)  found  that  the  interfacial  bond  strengths  and 
frictional  shear  strength  of  unfinished  PRD  with  nylon  12  was  6640  psi  and  292  psl, 
respectively. 

In  samples  fractured  by  cyclic  fatigue  (Figures  22  and  23^  the  aspect  ratio 
Is  lower  than  In  the  tensile  test  specimens.  Moreover,  one  can  see  In  Figure  23 
that  the  nylon  matrix  exhibited  a high  degree  of  yielding  and  a high  degree  of 
adhesion  to  the  fibers. 

In  the  PRD  reinforced  nylon  composite,  the  early  damage  occurs  near  fibers 
lying  perpendicular  to  the  stress  direction.  The  damage  took  place  in  the  form 
of  debonding  (l.e.  separation  between  the  fiber  and  the  matrix)  that  Intensified 
progressively,  affecting  fibers  at  smeller  angles  to  the  applied  load  until  at 
some  higher  load  (less  than  the  ultimate  load)  resin  cracks  occurred  In  the  resin 
rich  zone  of  the  material. 

For  graphite  reinforced  nylon  composites  a similar  mechanism  Is  suggested. 

The  </n!y  difference  Is  that  because  of  the  good  adhesion  between  the  matrix  and 
the  fibers,  the  resin  tends  to  draw  In  the  direction  of  the  stress. 
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VII.  CONCLUSIONS 

1.  Both  the  fatigue  lives  and  tensile  strengths  of  short  graphite  fiber 
reinforced  nylon  66  can  be  characterized  by  extreme  value  statistics.  Modified 
Wei bull  functions  may  be  used  to  describe  the  cumulative  distributions  for 
these  two  properties.  The  quality  of  a particular  fabrication  process  can  be 
evaluated  by  measuring  the  various  moments  of  these  distributions. 

2.  The  mechanism  of  fatigue  failure  in  graphite  fiber  reinforced  nylon  is  shown 

to  be  influenced  by  the  technique  of  fabrication  of  the  composite.  The  compression 
molded  materials  were  linearly  elastic  and  exhibited  an  isothermal,  brittle 
fatigue  failure  through  initiation  and  propagation  of  cracks.  The  injection 
molded  materials  were  non-linear  and  exhibited  considerable  ductile  flow  prior  to 
ultimate  failure.  Continuous  heating  of  the  material  resulted  in  premature 
failure,  suggesting  a thermal  mode  of  failure. 

3.  For  a specific  population  produced  by  one  fabrication  process  using  the 

same  raw  materials  for  each  member  of  the  population,  there  exists  a unique  relation 
between  the  cumulative  distribution  of  tensile  strengths  and  the  cumulative 
distribution  of  fatigue  lifetimes.  Thus  If  the  effect  of  an  environmental  variable 
on  the  distribution  of  strength  Is  measured,  the  effect  on  the  fatigue  lifetime 
can  be  estimated. 

4.  Since  the  rate  of  crack  propagation  In  graphite  reinforced  nylon  Is  highly 
sensitive  to  stress  (proportional  to  the  15th  - 20th  power  of  stress),  the 
prediction  of  fatigue  lifetime  directly  from  crack  propagation  data  Is  highly 
unrel  able. 

5.  Tne  mechanism  of  crack  propagation  In  short  graphite  fiber  reinforced  nylon 

Is  phenomenologically  similar  to  the  growth  and  fracture  of  crazes  In  thermoplastics. 
The  propagation  of  damage  at  a crack  tip  Is  controlled  by  matrix  deformation, 
cavitation,  fiber  breakage  and  fiber  pullout.  Damage  can  propagate  in  the  absence 
of  crack  growth  until  a critical  point  Is  reached  at  which  time  the  material 
fractures  catastrophically. 
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6.  Hybrid  composites  of  polyaramid  fibers  and  graphite  fibers  are  more 
resistant  to  crack  propagation  than  graphite  composites  with  an  equivalent  fiber 
content.  The  toughness  attributed  to  the  polyaramid  fibers  Is  translated  directly 
into  the  composite.  Increased  ductility  increases  the  sensitivity  of  the  fatigue 
crack  propagation  rate  to  the  amplitude  of  the  stress  intensity  factor  while 
decreasing  Its  sensitivity  to  the  maximum  stress  Intensity  factor. 
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During  Fatigue 
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Fig.  16.  Schematic  Structure  of  Hybrid  lamina  composite 
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Fig.  1$.  Crack  Propagation  Velocity  va  Maxima 
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TABLE  A-l 


Breaking  Strength  Distribution 
Population  Annealed  for  1/2  Hr 
Conditioned  at  55%  RH 


Fraction  Surviving  Breaking  Stress  (psi) 


t 

\ 


0.97 

9,450 

0.94 

10,242 

0.91 

10,247 

0.88 

10,868 

0.84 

10,868 

0.81 

11,000 

0.78 

11,173 

0.75 

11,401 

0.72 

11,750 

0.69 

12>  0C0 

0.66 

12,100 

0.63 

12,215 

0.59 

12,215 

0.56 

12,260 

0.53 

12,333 

0.50 

12,414 

0.47 

12,500 

0.44 

12,500 

0.41 

12,571 

0.38 

13,012 

0.34 

13,100 

0.31 

13,406 

0.28 

13,406 

0.25 

13,810 

0.22 

13,899 

0.19 

13,900 

0.16 

13,900 

0.13 

14,201 

0.09 

14,206 

0.06 

14,770 

0.03 

14,821 
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Table  a- 2 


Breaking  Strength  Distribution.Population  Annealed  for  1 Hour 

Conditioned  at  55%  RH 


Fraction  Surviving 

Breaking  Str 

0.94 

9465 

0.89 

10052 

0.83 

10572 

0.78 

10671 

0.  72 

10704 

0.67 

11084 

0.61 

11350 

0.56 

11669 

0.50 

12176 

0.44 

12413 

0.39 

12566 

0.33 

12588 

0.28 

12727 

0.22 

12725 

0.17 

12801 

0.11 

13220 

0.06 

14053 

Table  A-3 


i 


Breaking  Strength  Distribution.  Population  Annealed  for  1 Hour 

Conditioned  at  100%  RH 


Fraction  Surviving  Breaking  Stress  (p.  s.  i. ) 


0.  95 

3897 

0.90 

9445 

0.85 

9830 

0.  80 

10027 

0. 75 

10476 

0.  70 

10607 

0.65 

10645 

0.60 

10642 

0.55 

10701 

0.50 

10755 

0.45 

11240 

0.40 

11519 

0.35 

11550 

0.30 

11607 

0.25 

11628 

0.20 

11769 

0.15 

12055 

0.10 

12278 

0.05 

12723 
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Table  a- 4 


Breaking  Strength  Distribution  Population  Annealed  for 

Conditioned  at  55 % RH 


Fraction  Surviving 

0.96 
0.91 
0. 87 
0.  83 
0.78 
0.74 
0.70 
0.65 
0.60 
0.57 
0.52 
0.47 
0. 43 
0.39 
0.35 
0.30 
0.26 
0.22 
0.17 
0.13 
0.09 
0.04 


Breaking  Stress 

14103 

14200 

15016 

15177 

15177 

15933 

15985 

16000 

16000 

16080 

16160 

16200 

16295 

16410 

16480 

16570 

16700 

16989 

17081 

17090 

18760 

18845 


24  Hours 


(p.  s.  i. ) 


Table  a- 5 


/ 

i 

i 


t 

Breaking  Strength  Distribution  Injection  Molded  Population  Annealed 
for  24  Hours  Conditioned  at  55%  RH 


TABLE  A- 6 


Breaking  Time  Distribution  - A 1/2  RH55 
Specimens  Annealed  for  Half  an  Hour 
Conditioned  at  55%  R.H. 


max  = 11250  psi 
Zero  time  breaks  = 4 


P(t  ) 

t (Seconds) 

0.872 

100 

0.846 

130 

0.821 

160 

0.795 

160 

0.769 

200 

0.744 

230 

0.718 

340 

0.692 

620 

0.667 

630 

0.641 

830 

0.615 

1280 

0.590 

1290 

0.564 

1650 

0.538 

1680 

0.513 

1830 

0.487 

2070 

0.462 

2590 

0.436 

2610 

0.410 

2690 

0.385 

2980 

0.359 

*3010 

0.333 

4400 

0.308 

4580 

0.282 

5340 

0.256 

6020 

0.231 

6410 

0.205 

7850 

0.179 

8540 

0.154 

10430 

0.120 

13400 

0.10) 

17130 

0.077 

17440 

0.001 

39560 

0.000 

4482 

Breaking  Time  Distribution  - A1RH100-1 
Specimens  Annealed  for  One  Hour 
Conditioned  at  100%  R.H. 


“ 9061  psi 

Zero  Time  Breaks  = 1 
Infinite  Time  Breaks  = 4 


ptt.) 


tg  (Seconds) 


TABLE  A- 8 


Breaking  Time  Distribution-AlRH100-2 
Specimens  Annealed  for  One  Hour 
Conditioned  at  100%  R.H. 

= 9279  P51 

Zero  time  breaks  = 1 
Infinite  time  breaks  = 1 


tg (Seconds) 


.909 

.864 

.818 

.773 

.727 

.682 

.636 

.591 

.545 

.500 

.455 

.409 

.364 

.318 

.273 

.227 

.182 

.136 

.091 


84 

192 

229 

282 

312 

792 

1033 

1376 

1974 

2283 

2350 

5689 

5991 

8942 

14232 

33048 

34697 

79209 

98479 


TABLE  A- 9 


BreaJdng  Time  Distribution  - A1RH1Q0-3 
Specimens  Annelaed  for  One  Hour 
Conditioned  at  100%  R.H. 

^ max  - 9849  psi 
Zero  time  breaks  = 6 


P<V 

tg  (Seconds) 

0.696 

35 

0.652 

49 

0.609 

91 

0.565 

102 

0.522 

104 

0.478 

152 

0.435 

217 

0.391 

265 

0.348 

268 

0.304 

269 

0.261 

304 

0.217 

555 

0.174 

790 

0.130 

1157 

0.087 

1348 

0.043 

3513 

r 
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TABLE  A-10 

Breaking  Time  Distribution  -A24RH55 
Specimens  Annealed  for  24  Hours 
Conditioned  at  55%  R.H. 


* 13566  1*1 
Zero  time  breaks  = 4 


tg  (Seconds) 

p<v 

tg(Secc 

0.881 

81 

0.381 

3979 

0.857 

C8 

0.357 

4784 

0.833 

93 

0.333 

5479 

0.810 

102 

0.310 

5619 

0.786 

201 

0.286 

8145 

0.762 

244 

0.262 

8907 

0.738 

247 

0.238 

9877 

0.714 

618 

0.214 

13383 

0.690 

641 

0.190 

13505 

0.667 

707 

0.617 

19414 

0.643 

737 

0.143 

21308* 

0.619 

865 

0.119 

21482 

0.595 

909 

0.095 

21948 

0.571 

1041 

0.071 

42792 

0.548 

1220 

0.048 

49533 

0.524 

1385 

0.024 

92176 

0.500 

1995 

0.476 

2581 

0.452 

2801 

0.429 

3042 

0.405 

3709 

TABLE  A- 12 

Data  for  Crack  Propagation  Rate  Plots 
A1  RH  100 


154-2 


154-3 


154-5 


154-1 


138 


| * iff* 


K 


max 


(1-R) 


0.6 


ln(&a/  At) 


1 

5.90 

-6.  78 

2 

5.  91 

-5.  99 

3 

5.  92 

-4.  95 

4 

5.  94 

-4. 42 

5 

6.03 

-6.91 

6 

6.  05 

-5.  96 

7 

6.  03 

-2.65 

8 

6. 22 

-2.  53 

9 

6.25 

-1.  63 

10 

6.27 

0.  6413 

11 

5.  87 

-5.  86 

12 

5.  90 

-6.78 

13 

5.  94 

-5.46 

14 

6.04 

-2.  93 

15 

6.  07 

-1.  00 

16 

6.23 

1.9/ 

17 

5.85 

-5.28 

18 

5.  96 

-5.37 

19 

6.  00 

-5.20 

20 

5.  99 

-3.91 

21 

6.  08 

-4.14 

22 

6.11 

-2.81 

24 

6.12 

-4. 72 

25 

6.14 

-3.17 

26 

6.19 

-2.  00 

27 

6. 22 

-1. 44 

28 

6.25 

0.06 

TABLE  A- 13 


I: 

-I 

[ \ Data  for  Crack  Propagation  Rate  Plots 

\ A1  RH  55 


v- 

max 

(1-R)0, 1 

In  (a  a^A  t) 

O 156-3 

1 

6. 11 

-6.  66 

‘ 

2 

6.20 

-6.  03 

3 

6.26 

-5.60 

i 

4 

6.29 

-3.  04 

< 

‘ 

5 

6.31 

-1.  85 

156-4 

6 

6.24 

-4.42 

* i 

7 

6.25 

-3.  65 

f 

8 

6. 31 

-5.  08 

< » 
* 

9 

6.33 

-3.25 

10 

6.37 

-2.  51 

■/  1 

11 

6.38 

-1.  76 

< 1 

12 

6.45 

-0.  419 

4 

156-2 

13 

6.  07 

-6.  68 

14 

6.16 

-6.  52 

15 

6.21 

-3.64 

i 

16 

6.25 

-4.61 

\ 

156-5 

17 

6.  06 

-5.  90 

I 

18 

6.14 

-6.80 

19 

6.20 

-5.77 

I 

20 

. 6.23 

-5.78 

f 

21 

6.26 

-4’.  54 

1 

22 

6.27 

-2.19 

} 

% 


5 


! 


t 

II 

TABLE  A- 14 


Data  for  Crack  Propagation  Rate  Plots 
A 2 4 RH  55 


TABLE  A- 15 


Apparent  Fracture  Toughness  Obtained  at 
Different  R Values  During  Crack  Propagation  Experiments 


Sample 

Experiment 

rr 

app 

R 

A24RI1  55 

127 

912.  2 

0.  42 

129 

988.2 

0.18 

143 

1077.6 

0.43 

144 

1263.2 

0.  73 

AIRH  55 

156-2 

576.2 

0.2 

156-3 

602.2 

0.4 

156-4 

670.  0 

0.11 

156-5 

606.  8 

0.  5 

634.4 


AIRH  100 


141 


0.14 


TABLE  A- 16 


Threshold  Values  Observed  During 
Crack  Propagation  Rate  Experiments 


Population 

Experiment 

V 

max 

R 

In  ERD 

In  V 

A1RH100 

154-2 

493.9 

0.419 

5.58 

-6.07 

154-3 

440.8 

0.146 

6.10 

-6.22 

A1RK55 

156-5 

429.8 

0.412 

6.01 

-6.65 

A24RH55 

126 

793.6 

0.52 

6.38 

-6.91 

144 

1018.2 

0.714 

6.43 

-8.11 

14  3 

915.1 

0.429 

6.59 

-9.21 

l 


Table  a- 17 


» 


Data  for  aB(L.)  vs  t (L.)  Curve 
A 1/2  RH55 


p 

log  tB 

log  <Jj 

ai 

4.176 

4.158 

0.2 

3.  929 

4.140 

0.3 

3.  724 

4.127 

0.4 

3.  580 

4.114 

0.5 

3.342 

4.102 

0.6 

3.114 

4.  088 

0.7 

2.  748 

4.  074 

0.8 

2.  255 

4.  055 

A1RH100 

• 

P 

l°g°B 

log  *B-1 

log  tB_;. 

!°g  tB_ 

0.05 

4.104 

3. 230 

0.1 

4.  092 

4.  643 

3.  079 

0.2 

4.  076 

4.  869 

4.260 

2.  845 

0.3 

4.  064 

4.301 

4.0 

2. 653 

0.4 

4.  057 

3.964 

3.732 

2. 447 

0.5 

4.  045 

3.653 

3.477 

2.190 

0.6 

4.033 

3.362 

3.161 

1.  833 

0.7 

4.021 

2.  954 

2.799 

0.8 

4. 003 

2.477 

2.362 

0.9 

3. 973 

1. 778 

1. 531 

TABLE  A- 18 


i 


Shifted  Data  for  ^ (L.)  vs  (7b  (L^)  Master  Curve 


log  <tB 

>s 

10BVo 

p 

log 

AIRH100-1 

AIR  HI  00-3 

AIRH100-2 

0.1 

4.  092 

4.233 

4.  643 

0.2 

4.  076 

4.  332 

3.904 

4.260 

0.3 

4.  064 

3.797 

3.647 

4.0 

0.4 

4.  057 

3.494 

3.375 

3.  732  ' 

0.5 

4.  045 

3.215 

3.  054 

3.477 

0.6 

4.  033 

2.  964 

2.  618 

3.161 

0.7 

4.021 

2.601 

2.799 

0.8 

4.  003 

2.179 

2.362 

0.9 

3.  973 

1.  585 

1.  531 
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TABLE  A- 19 


Generation  of  a Single  Fatigue  Distribution 
Curve  (Fig.  29)  by  the  Use  of  Crack  Propagation 
Law  and  Breaking  Strength  Distribution 

Starting  Point:  ^ max  = 11250  at  ERD  = ERD  crit.  = 635 
• P ( ^ max)  = 0.815 


p 

max  PSI 

ERD 

•cm 

a— i — 
nun 

tD  (sec) 

0.1 

14,400 

496.1 

3.14xl0“3 

8 . 2x1.  (P 

0.2 

13,800 

517.7 

4.73xl0-3 

5 . 4xl03 

0,3 

13,400 

533.1 

6.36x10-3 

4.05x103 

0.4 

13,030 

548.3 

8.62x10-3 

3.0x103 

0.5 

12,650 

564.7 

1.227x10-2 

2.1x103 

0.6 

12,290 

581 . 3 

1.84x10-2 

1.4xl03 

0.7 

11,860 

602.3 

3.55x10-2 

7.26x102 

910773  2.149372  134.9354 


45.0  12.5  57.5  8.171281  3.552731  53.15021 


40.0  15.0  55.0  2.373295  12.945524  15.31301 


0.805  50.0  10.0  60.0  6.824553  2.274851  55.40770 


FATIGUE  CRACK  PROPAGATION  RATE  DATA  FOR 
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